Landau singularities and singularities of holonomic 
integrals of the Ising class 



S. Boukraa^ S. Hassani § , J.-M. Maillard* and N. Zenine § 

fLPTHIRM and Departement d'Aeronautique, Universite de Blida, Blida, 
Algeria 

§Centre de Recherche Nucleaire d'Alger, 

2 Bd. Prantz Fanon, BP 399, 16000 Alger, Algeria 

J LPTMC, Universite de Paris 6, Tour 24, 4eme etage, case 121, 

4 Place Jussieu, 75252 Paris Cedex 05, France 

E-mail: maillard91ptmc.jussieu.fr, maillardSlptl.jussieu.fr, 
njzenine9yahoo . com, boukraaOmail .univ-blida.dz 

Abstract. 

We consider families of multiple and simple integrals of the "Ising class" and 
the linear ordinary differential equations with polynomial coefficients they are 
solutions of. We compare the full set of singularities given by the roots of the 
head polynomial of these linear ODE's and the subset of singularities occurring 
in the integrals, with the singularities obtained from the Landau conditions. For 
these Ising class integrals, we show that the Landau conditions can be worked out, 
either to give the singularities of the corresponding linear differential equation or 
the singularities occurring in the integral. The singular behavior of these integrals 
is obtained in the self-dual variable w = s/2/(l + s 2 ), with s = sinh(2X), where 
K = J/kT is the usual Ising model coupling constant. Switching to the variable 
s, we show that the singularities of the analytic continuation of series expansions 
of these integrals actually break the Kramers- Wannier duality. We revisit the 
singular behavior (J. Phys. A 38 (2005) 9439-9474) of the third contribution to 
the magnetic susceptibility of Ising model \^ a * the points l + 3to + 4u> 2 = and 
show that x' 3 ^ («) is not singular at the corresponding points inside the unit circle 
\s\ = 1, while its analytical continuation in the variable s is actually singular at 
the corresponding points 2 + s + s 2 = oustside the unit circle > 1). 
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1. Introduction 

Understanding the magnetic susceptibility \ of the Ising model as a function of the 
temperature remains an outstanding problem in lattice statistical physics. One of 
the problems raised in this model is the possibility of occurrence for x of a natural 
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boundary in the complex plane of the variable s = smh(2K), where K = J/kT is 
the conventional Ising model coupling constant. There is a strong argument on this 
point, since it has been shown by Nickel [2HI that the singularities of each n— particle 
contribution x^ to the susceptibility are given by 

s + - = u k + u- k +u' n + u- m , (1) 
s 

with: u n = 1, — n <k, m < n 

which are actually lying on the unit circle \s\ = 1. In the following , we will call 
these singularities "nickcllian singularities" . These individual x are n ~ particle 
contributions to the infinite sum [3] 

X = £ X {n) (2) 

where the singularities |T]) accumulate and (in absence of cancellation) density the 
unit circle, possibly creating a natural boundary for x- The above sum is restricted 
to n even (respectively odd) for the low (respectively high) temperature. 

Among these n— particle contributions, x^ an d X^ were known long ago, and 
the following two, x^ an d X^\ have been shown recently [H[5j [6], to satisfy Fuchsian 
linear differential equations of, respectively, order seven and order ten. 

Our study is based on a holonomic approach which uses the fact that, the 
integrand being an algebraic function of all the variables, the multiple integral under 
consideration, can be shown to be holonomic [7], i.e. it satisfies finite order linear 
differential equations, with polynomial coefficients, in the remaining variable for which 
no integration has been performed. The analysis, then, focuses on the corresponding 
linear differential operators. This holonomic approach gives the full set of singularities 
as roots of the head polynomial in front of the highest derivative of the differential 
equation (excluding the "apparent" polynomial [HOIS]). 

The linear ODE corresponding to x^ nas 0, besides the known singularities 
s 2 = 1, and s 2 = —1, only the singularities |T]) predicted by Nickel [HH], while the 
linear ODE corresponding to x^ 3 ) has shown a pair of singularities 1 + 3w + Aw 2 = 
(where w — s/2/(l + s 2 )) that are not on the unit circle \s\ — 1. The condition \s\ — 1 
is equivalent, for the variable w, to be real and such that its absolute value is greater 
than 1/4. In this paper, by w is "in W c ", we will mean w e] — oo, —1/4] U [1/4, oof. 

Therefore, the following question is in order. Should we expect in the linear 
ODE's corresponding to the successive n-particle contribution x^ n \ the occurrence of 
singularities that are of similar nature as 1 + 3iu + Aw 2 — 0, i.e. not lying in W c , or 
only the occurrence of the "nickellian singularities" given by |T])? 

These expected singularities of the linear ODE governing the x^ are n °t 
necessarily singularities of the particular solution x ■ Do singularities with |s| > 1 
occur in the X i an d accumulate, the question of a natural boundary for x becomes 
of the most importance. 

In previous works, we have introduced a "connection matrix method" [S] to 
determine if a singularity of the linear ODE is actually a singularity occurring in 
the multiple integral. For the three-particle contribution x^ 3 \ we have shown [8] 
that the two quadratic roots of 1 + 3w + Aw 2 = 0, singularities of the linear ODE, 
are actually absent in the multiple integral defining x^ written as a function of the 
self-dual variable 

J Be aware that some subtlety can be encountered when the variable s is used, see Section (3.3) and 
Section (4.5). 
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In our holonomic approach, the singularities of the multiple integrals are, thus, 
obtained in two steps: find the linear differential equation, then use our "connection 
method" 8 in order to analyze the singularities of the solutions of this ODE and, 
especially, the singularities of that particular solution corresponding to the multiple 
integral. These computations may be rather cumbersome, depending on the order of 
the linear ODE, on the number of singularities, and their distribution in the complex 
plane of the variable. 

In particle physics, in the study of analytical properties of the S-matrix [9j, or in 
Feynman diagram calculations, the amplitudes of the Green functions are defined 
through an integral representation and depend on external kinematical invariants 
(the Mandelstam variables), the quadri-momentum in the loops being the integration 
variables. The singularities of these integrals are associated with physical thresholds 
in each of the physical regions in which they describe the process under consideration. 
The question of forecasting the location of the singularities of functions defined through 
integral representations is, thus, a general and interesting problem in mathematical 
physics [TO] . 

An analytical approach [9] of this issue provides, by imposing conditions on 
the integrand, singled-out sets of complex points, as possible singularities of the 
integral. These conditions on the integrand are called Landau conditions [9] 111] 1 1 2 j . 
The singularities obtained this way are called Landau singularities. Note that, 
the Landau singularities were introduced to be those of the integral, and that, 
the Landau conditions are only necessary conditions for the existence of such 
singularities [9l[l2l[T3]. In fact, when this analytical ^-matrix approach was flourishing 
in the litterature, it was not yet known [7j that the integrals of (physical) interest were 
actually holonomic, that is, solutions of finite order linear differential equations with 
polynomial coefficients. The fact that the Landau singularities could have a better 
overlap with the singularities of the linear ODE's, rather than with the singularities 
of the integral of interest, was a question nobody could imagine at that time. 

In many domains of mathematical physics [14] . the computation of multiple 
integrals like the evaluation of Feynman diagrams 10J, or various correlation functions 
in quantum field theory and statistical mechanics [T^l UJ)] is an important problem. 
The n-particle contribution to the magnetic susceptibility of the Ising model have 
integral representations. They are given by (n — l)-dimensional integrals (2] [17\ ITS] 
that read (omitting the prefactor) 

1 n— 1 /.27T ii n o 

h-(Ul £)(n«0 ■*">■(*"')■ < 3 > 



3=1"" 3=1 



where 



and 



witr§ 



G{n)= (n^) ( "- i)/2 . n 2sin | ( ^ )/2) > w 

Ri n) = i + feu xi i (5) 



i-n 



2w 



i) + ] /{l-2wcos(&)f ~4w 2 



(6) 



1 — 2w cos( 



% Variables Xi and yi are the variables noted Xi, and yi, in previous papers [3][5][6]. 
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2w 



■n, 



Vi = 



(1 - 2uicos(0i)) 2 - 4w 2 



(7) 



valid for small w and, elsewhere, by analytical continuation. Recall that small w 
corresponds to small values of s, as well as large values of s. 

If successful, the Landau conditions will give the singularities of these integrals 
bypassing the linear ODE search and the singularity analyzis from our connection 
method [H], possibly providing a deeper insight on the natural boundary problem of 
the susceptibility. 

Unfortunatly, the multiple integrals ([3]) are too involved for the Landau conditions 
to be worked out. We will consider, instead, and as a preliminary project, two 
kinds of functions in the form of integral representations (one-dimensional and 
multidimensional) which belong to the Ising clasi^J. The integrands of this class of 
integrals have a trigonometric-hyperbolic form and are simple algebraic functions in 
the cosinus of the angular integration variables, and are also algebraic in the remaining 
w variable. These integrals are thus holonomic: they must satisfy finite order linear 
ODE's with polynomial coefficients in the remaining variable w. The route we take 
for obtaining the linear ODE's of the two sets of integrals will be through series 
expansions [HOIS]. 

One example corresponds to a multiple integral, but yields quite simple sets of 
singularities, whereas the second example corresponds to simple integrals, but yields 
quite non trivial sets of singularities. These two "Ising class" examples have been 
chosen to provide an "equivalent level of complexity" but complementary results on our 
singularity problem. For each example, we compare the singularities of the linear ODE, 
the singularities of the integral, and the singularities given by the Landau conditions. 

The paper is organized as follows. In Section ©, we recall the Landau equations 
method [9l [19] . Section ([3]) presents the first "toy model" for which we can obtain for 
any value n of this n-multiple integral, the corresponding linear ODE, remarkably 
of second order, and the corresponding solutions. Section (j4|) introduces a second 
"toy model" of simple integrals, gives the corresponding linear ODE's with their 
singularities, presents the Landau singularities, and compares these two sets with the 
actual singularities of the integral. In both sections, we consider some subtleties that 
appear when we switch from the self-dual variable w to the variable s, with respect 
to the analytic continuations of each. We finally comment the calculations on Landau 
conditions, sketch the forthcoming "Ising class" integral calculations we will study, 
and conclude. 

2. Landau conditions 

For the paper to be self contained, we recall, in this section, general features about 
the Landau singularities. Let us, for instance, consider the following integral: 



where f(w, <f>\, • ■ ■ , <f> n ) denotes an algebraic function of w and of trigonometric 
functions of the angles fa. Denote by Bj (w, <j>), j — 1, 2, • • • d, the set representing the 
boundary of the integration domain, and by Sj (w, cf>), j = 1, 2, • • • s, the set of varieties, 




(8) 



jj The terminology integral of the Ising class has also been used by Bailey, Borwein and Crandall in 
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locus of the singularities of the integrand (<p denotes <pi, <j>%, ■ • • , Landau equations 
amount to finding the necessary conditions for the singularities to occur when the 
hypercontour is pinched between the surfaces of singularities (pinch singularity) or 
meets a boundary variety (end-point singularity) . For this to happen, the parameters 
at and Pi should exist, not all equal to zero, and such that at the point (w, 4>) 

cxj ■ Bj{w,(t)) =0, j = 1, d (9) 
(3j ■ Sj(w, <t>) = 0, j = 1, • • • , s 

'Sr(^2 ai ' Bi ( w -> 4>)+^Pi- Si(w,(f>)j = 0, j = 1, • ■ • , n 

^ i i 

Solving this set of equations will give the singularities w that are "candidates" to be 
those of the integral ([8]) . Note that the easy resolution of these equations will depend 
on the number of integration variables and the number of singularity varieties. 
Consider, for instance, a one-dimensional integral in the form: 

" f D^l < 10 > 
Jo D(6, w) 

where D(8, w) represents the variety of singularities. From @, the Landau conditions 
readjj: 

a- D{6, w) = 0, (3-9=0, 

^(a-D(6,w) +0-O)= (11) 

Knowing that a and (3 should not be both equal to zero, it is easy to see that the set 
of equations (jTTJ) gives two kinds of solutions, namely: 

D(6, w) = 0, 6 = (mod 2tt) (12) 

or: D{0, w) = 0, ]fi D ( > w ) = ( 13 ) 

The first solution corresponds to endpoint singularities while the second one 
corresponds to pinched singularities. Note that mixed situations with endpoint 
together with pinched singularities, can occur for multidimensional integrals. 



3. A multiple integral: the y 2 product 

Recall that a multidimensional integral of the product of the yt quantities has been 
recently introduced by Nickel [13] to suggest that the singularities, like 1 +3w +4w 2 = 
0, appearing in [4] for the linear differential equation are pinch singularities. 

In this section, for elegance, and less prohibitive calculations purposes, we 
consider, instead, a multidimensional integral of the product of yf''s. The results 
are similar to the product of the y^s, as far as the singularities location is concerned 
(see Appendix A for the linear ODE bearing on the product of y\ with n = 3). 

Our multidimensional integral thus reads: 




% Strictly speaking, one has to add j3{6 — 2ir) = 0, but this equation does not change the discussion. 
We are considering 2n periodic integrands. 
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with 

n 

X>=0 (15) 

k=l 

To find the linear differential equations for these integrals, we make use of the 
Fourier expansion of yf which reads (dropping the indices) 

2 00 

£2 = £>(Q) + 2 D ( k ) ■ cos ( fc ^) 
fc=i 

00 00 
= £ D(k)Z k = £ £(fc).Z- fe (16) 
fc— — 00 fc— — 00 

with Z = exp(i <j)) and 

D(k) = D(-k) = .d(|Jfc|) (17) 

where d{k) is a non terminating hypergeometric series given by 

d(k) = (1 + fc) x (18) 
^ /3 fc k 5 k 3 fc3 , 3 , „„ « 
4F3 U + 2' 1+ 2'4 + 2'2 +2 ; 2' 1 + fc '2 +fc;16w ' 
The Fourier expansion (fT6)) . together with the definitions (fT7| . lfl8]l. allow to write the 
expansion of the integral Y^(w), fully integrated over the angles, as (see Appendix 
B for details) 

00 

F (n) (w) = (4w 2 )" ■Y,C k -w nk -d(k) n (19) 

fe=0 

where Cq = 1 and Ck — 2 for fc > 1. 



5. i. Linear ODE's and their solutions 

Series generation using (fT9|) can be carried out without any difficulty. We have 
computed the linear differential equations satisfied by the yW's up to n = 10. These 
Fuchsian linear ODE's have many remarkable properties. They are of order two, 
independently of the dimension n of the integral. All these Fuchsian linear ODE's 
have the singularities w = 0, w = 1/4, w = — 1/4 and 00, with respectively the critical 
exponents (2n, n - (-1)™), (-(n- 2)/2,-(n+ l)/2), (-(n - 2)/2, -(n + (-l) n )/2) 
and (—(71— l)/2, — (n + l)/2). The singularity w — is an apparent one. These linear 
ODE's have other regular singularities with critical exponents —1 or 0. 

Let us give, as an example, the linear ODE for n — 3 {F{w) denotes Y^\w)) 

2 d m 

E a rn(w)-—F(w) = (20) 
m=0 

where 

a 2 (w) = w 2 -(l-w) (l + 3w)(l-16w 2 ) 2 

x (1 - 5 w) (l + 3 w + 4w 2 ) • P 2 (w), 
ai(w) = (l — 16uj 2 ) w ■ Pi(w), ao(w) — 6Pq(w) 



Landau singularities of holonomic integrals of the Ising class 



7 



the polynomials being: 

P-2 = 1 - 2 w - 24 w 2 - 5 w 3 + 19 w 4 + 98 w 5 - 574 w 6 - 2835 w 7 
+ 2330w 8 + 392u; 9 , 

Pt = -9 + 16 w + 410 w 2 - 8 w 3 - 5164 w 4 - 7196 w 5 + 12896 w 6 
+ 65824 w 7 - 192637 w 8 - 1302860 w 9 + 327886 w 10 
+ 6666128 w 11 + 3469546 w 12 - 6576560 w 13 



2874016 u; 14 + 1967104-, 



,15 i 7COC4n „,,16 



P = 4 - 6 w - 224 w 2 + 29 w 3 + 3808 w 4 + 3272 w 5 - 21680 w 6 

- 71496 w 7 + 87964 w 8 + 1 146670 w 9 - 856344 w 10 

- 9773745 w 11 + 2595240 w 12 + 30066076 w 13 + 1626384 w 14 

- 26569152 w 15 - 8783616 w 16 ~ 1204224 w 17 + 2007040 w ls . 

The roots of polynomial P2(w) are apparent singularities. The linear ODE for n = 4 
is given in Appendix C. 

This second-order linear differential equation can be solved exactly, and has the 
general solution (ci and C2 are arbitrary constants): 

F(w) = c 1 -F 1 {w) +c 2 -F 2 {w), (21) 

Fi(w) = 

(3 - 59 w 2 - 82 w 3 + 386 w 4 + 8 w 5 + 64 w 6 ) • w 4 
3 (1 - Awf (1 - w) (1 + 3 w) (1 + Aw) (1 - 5 w) (1 + 3 w + 4 w 2 ) ' 
(l + 2w - 2w 2 ) • w 4 

p 2 H - 



(1 - 5 iu) (1 + 3 w + 4 w 2 ) VI - 16 

Each component Fi(w) is a solution of an order-one linear differential equation. All 
the singularities are poles, except w = ±1/4 which are branch points for the second 
component F2(w). For arbitrary constants of integration, the general solution F(w) 
has the same singular behavior. To obtain the particular solution corresponding to 
the integral (|14p for n = 3, i.e. F(w) — Y^ 3 \w), the constants of integration have 
to be fixed to c\ = 12, c 2 = —c\. Note that, when these constants of integration are 
fixed at these values, the solution F(w) — Y^ 3 \w) is no longer singular at w = 1/5, 
and at the quadratic roots wq of 1 + 3w + Aw 2 = 0, as can be easily checked: 

37 

Y H= TI^o"' when (22) 

v(3)r ^ 4099 , .3881^7 

F W= 90112 ±i ^0Tl2-' Wh6n W ^ W0 

However, these polynomials l-5w and 1 + 3 w + Aw 2 still remain in the expression 
of the solution. Note that the corresponding singularities 1 - 5w = and 
1 + 3w + Aw 2 = are not in W c , while the others in (|21|) are. 

As a second example, consider, for instance, the case n — 5. The singularities of 
the linear ODE, besides 1 — 16w 2 — 0, correspond to the roots of the two polynomials: 

Qf= (l-3w + w 2 ) (l + 5w + 5w 2 ), (23) 
Q { 2 5) = (l + 5u; + 13w 2 ) (l-7w + 5u; 2 -4u; 3 ) 

fl - 8 w + 20 w 2 - 17 w 3 ) (1 + 8 w + 20 w 2 + 15 w 3 + A w 4 ) . 
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One can check that the set of values of w such that 1 — 16w 2 = 0, and the singularities 
of the set = 0, are indeed singularities of the multiple integral Y 1 - 5 ^ (see below 
and Appendix C), while those of the set — are not. Here also, the singularities 
cannot be "simplified" due to the presence of the radical yl — 16u> 2 . One can also 

(5) 

check that the singularities corresponding to Q\ — are in W c , while those of the 
set Q 2 5) = are not. These remarkable facts hold for all the other Y^'s. 

With the differential equations up to n — 10, we find that the singularities are 
roots of the following polynomials given by Nickel [13] 

cos(fc0+) - cos ((ra - fc)<£ - ) = 0, (24) 

cos (0+) = 1 + i_, cos (J>- ) = -1 + i (25) 

which can be written as 5J: 

fj_ ^ „ /I 



Tk VM, + *J ~ Tn - k VhL ~ 1 ) = °' fe = 0, 1, ---.n (26) 
where Tk(w) are the Chebyshev polynomials of first kind. Note that, in this model, the 
singularity 1 + Aw = is not given by this formula for n odd. From these singularities, 
only those corresponding to k — 0, and k = n, are singularities of the Y^'s, i.e. 

1 - T »(i- 1 ))-( 1 - T "(i + 1 ))= °- < 27 > 

We are, now, in position to write the generic structure of the solutions F' n '(w). 
For n odd, the solutions 7'"'(ai) have the form: 

Y (n \w) = (28) 

c i" } ■ M n) _l» d(») 



(1 - 16 w 2 ) 11 ' 2 - 1 ■ Q { 2 n) Vl-16w 2 (l - 4w) • 
For n even, they read: 

n— 1 

y (n) H 



c 2 ' ^2 



((1-16U. 2 )™/ 2 - 1 ) -Q| n) H 

4 n) .flWw 4 n) -fl (n) (-«0 



V(l +4w) 3/2 • Qi" } H (1 - 4w) 3/2 • 0[ n) (-w) / 

The coefficients and c^""* are integration constants. The polynomials R^ n \ 
and i?2™^ > m these expressions, have rational coefficients and are such that the constant 
coefficient in w is equal to 1. The polynomials Qi and Q 2 "' are constructed, 
respectively, from the roots of (f27|) and from the roots of the cumulative product (for 
k = 1, • • • , n — 1) of ([26]) . The set of roots = (resp. Q2 = 0) are (resp. are 

not) in W c - 

The general solutions above give the particular solution (fT4| for — — . 
When the integration constants satisfy this condition the solution is no longer singular 
at the roots of . 

% Note that Nickel 1131 excludes the cases k = and k = n. 
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3.2. Landau conditions for 

Consider, now, the Landau conditions^ for these integrals which amounts to 

solving 




for otj and f3j not all equal to zero. The Landau singularities are obtained for 
the following configurations. In the first configuration, a n = 0, one obtains the 
singularities w — ±1/4. The second configuration ay ^ 0, for all j, gives when at 
least one j3j is not zero, the singularity w = 1/4. For j3j — 0, for all j, one obtains the 
equations ([25]) for, respectively, k angles <fi + and n — k angles <f>~ . The constraint (|15p 
becomes: 



Using only the first condition (j3"Tj) . one obtains (|26p which are the singularities of 
the linear ODE, while using both conditions (|31[) . (j3"2"|) , one obtains the polynomials 
corresponding to k = and k — n in (|26[) . which only give the singularities of the 
multiple integral. The last configuration is embedded in the second configuration 
when aj = for j = 1, 2, • • • ,p and aj ^ for j = p-\- 1, • • • , n. These polynomials are 
just given by (130p with n~p instead of n, since, in the constraint (|15|) . the angles (^j, 
j = 1, • • • ,p become equal to zero by the equations (|29j) . The Landau conditions for 
give then the singularities of Y^ p \ with p — 1,2, • • • ,n — 1. These singularities 
are spurious with respect to the linear ODE. 

5.5. Kramers- Wannier duality subtleties: w versus s 

In the integrals of the previous section, the calculations are performed in the variable 
w. In particular we have made no distinction between the high and low temperature 
regime (s small or large). When trying to bridge the families of "Ising class" integrals 
we are presenting in this paper and the actual Ising model integrals x*"\ we should 
point out the subtleties of using the self-dual variable w which does not distinguish 
between the interior of the unit circle \s\ < 1 and the exterior of the unit circle |s| > 1. 

Let us come back to solution , and let us consider the singularity 1 — 5w = 0. 
The solution Y^iw) is not singular at w = 1/5 where it takes a finite value (|22p , 
although the polynomial 1 — 5w still remains in the algebraic expression of the function. 
One has a similar situation for the singularities 1 + 3w + 4w 2 = 0. The non 
"cancellation" of these polynomials, in the algebraic expression of the solution, may 



k <p + + (n 
which is equivalent to 



k)(f>- = 2tt 



(30) 




(31) 
(32) 



H See the resolution of Landau conditions for the more elaborate family of integrals of Section (4). 
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mean that the term 1 — 16u> 2 , occurring in the square root (|21|) , is a perfect square 
when written in some other variable. The last is nothing else than the variable s. 

The analytic continuation in the variable s is performed by writing the linear 
ODE corresponding to in the variable s, and solving. One obtains the following 
solutions 

F(s) = m-Y^is) + n 2 -Y^(^j where (33) 

Y^(s) = ^ with: 

N Y = (6s 7 + s 5 -lls 4 + 8s 3 -32s 2 + 16s-16) • s 6 , 
D Y = (l + s 2 ) 2 (s + l) 2 (s-l) 4 (2s 2 + 3s + 2) (2 s 2 - s + 2) 
x (s 2 + s + 2) (s - 2) 

The polynomials (1 — 5w) and (1 + 3w + 4w 2 ), written in the variable s, 
factorize, respectively, as (1 — 2s) • (s — 2) and (1 + s + 2s 2 )(2 + s + s 2 ) . The 
"physical" solution y( 3 )(s) is not singular at the point 1 — 2s = and at the 
quadratic points 1 + s + 2s = 0, which are inside the unit circle (|s| < 1). Do 
note, however, that the "physical" solution Y^(s) actually presents a singularity at 
s — 2 = 0, as well as at 2 + s + s 2 = 0, which are outside the unit circle (|s| > 1). 
The same conclusions hold for FW, n > 3. 

Let us now present a second family of more elaborate integrals for which we 
address the same questions: do the Landau singularities identify with the singularities 
of the linear ODE or with the singularities of the integral ? Do the subtleties 
encoutered using the self-dual variable w instead of the s variable also occur in this 
family? 

4. The Diagonal model 

Let us first recall the fundamental, and deep, mathematical relation between 
Hadamard product of series, multiple integration, and diagonal reduction |20[ I21j . 
Christol's conjecture amounts to saying that any algebraic power series in n variables 
is the "diagonal" of a rational power series of 2n variables ;2Q, 22, 23, 24,. These kind 
of quite remarkable results [24[ [25l [26l [27] gave us the idea to build our second "toy 
model" by considering the "diagonal" of the integrand (seen as a function of (n — 1) 
angles <j>j) occurring in the multiple integrals corresponding to the X s i an d thus 
perform, instead of n — 1 integrals, only one integration. 

At first sight, one can imagine that such a drastic reduction completely 
trivializes the function one tries to calculate, and yields a new function with totally 
different singularities. On the other hand, in view of the previous non-trivial 
results 20, 22} [23l EU [25l [26l [27], one can, a contrario, imagine that such a "diagonal 
reduction" procedure keeps "some" relevant analytical informations, especially on the 
singularities. Replacing the quite involved multiple integral by an integration on only 
one angle, and since the integrand is an algebraic function of trigonometric functions 
of this angle and of the variable w — s/(l + s 2 )/2, the result will be holonomic. 
This integral on only one variable will be solution of a finite order linear differential 
equation in w with polynomial coefficients. Could it be possible that the singularities 
of this (probably Fuchsian) linear ODE, keep some "memory" of the singularities of 
the involved initial multiple integral? 
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The second set of integrals, we present in this paper, make use of the following 
simplifications. From the integrand of m © 5 quicking out and the product 
on we only keep the following quantity in the integrand: 

i.i+nkf* = i +i i 

n\ i-nr=i^ « ! - ! i-nr=i^ 1 j 

This is reminiscent of the assumptions made in pQ that the singularities (or a subset 
of them) arise from the symmetry points of the integrand and the vanishing of the 
denominator of the above quantity. 

As previously explained, instead of performing the integration on (n — 1) angles 
(f>j, we integrate the above quantity on the principal diagonal, i.e. from a (n — 1)- 
dimensional integral we reduce ouselves to only one integration angle <p : 

4>x = (f>2 = ■ ■ ■ = 4>n-i — 4>, 4>n — -(n - 1) 4> (35) 

Call the integral on cf> of $M$: 

-(») = _1 lf 2 ^ 1 fo^ 

D n\ n\ J 2tt 1 -x"- 1 ^) ■x((n- 1)0) ( ' 

where x{<p) is given by ((6|). Expression (|36|) . when expanded in the variable x, 
becomes: 

The term x p has a Fourier expansion given by [5]: 

OO 

x p = w p ■ (b(0, P ) +2^2 w k -cos kcj)-b(k,p)^ (38) 
fc=i 

where b{k,p) is a non terminating hypergeometric series that reads (with m = k + p): 

b(k,p) = ("^ x (39) 
' (1 + m) (1 + m) (2 + m) m 



t*3( 2 2 V 2 ; ,y;l + fc,l+p,l + m; 16w 2 J 



With this Fourier series, the integration rules are straigthforward. 

T 
1 



(n) 

The fully integrated expansion of $> D simply reads: 



= -- (40) 



oc oc 



- EE W n[k+P) ■C k -b{k,p)-b(k{n~l),p{n-l)) 



p=0 k=0 



where Cq = 1 and Ck — 2 for fc > 1. These sums are clearly less consuming in the 
computational effort as compared to the multidimensional integrals ([3]). 



(n) 

4-1. Fuchsian linear differential equations for <f> D and their singularities 

(n) 

Generating sufficiently large series in w for each <£ < D , we have actually been able to 
find the corresponding linear differential equations up to n — 8. The linear ODE's for 
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, are of order four and are given in Appendix D. For instance, the linear 
ODE satisfied by reads: 

4 d m 

E a -N-^r F H = ( 41 ) 

m— 1 

where 

a 4 (w) = (1 - w) (1 + 2w) (l - 16ui 2 ) (1 +3w + 4w 2 ) w 2 ■ P 4 (w) 

a 3 (w) = w ■ P 3 (w), a 2 (w) = P 2 (w), ai(w) = Pi(w) 

the polynomials Pi(w) being displayed in Appendix D. 

One does remark, on this second example, that the singularities of the linear ODE 
for IS IS] actually pop out in the linear ODE of this "toy integral" , in particular 
the quite non trivial 1 + 3u> + 4u> 2 = singularities. 

The linear differential equation (|4"Tj) is of Fuchsian type. It has just regular 

singularities: one can actually see that the degree eleven polynomial P4 does 

(3) 

correspond to apparent singularities. The linear ODE's for $> y D ' and for the next 

s, all have logarithmic singularities at the points w = 0, ±1/4 and w = 00. The 
other regular singular points have a x~ x l 2 singularity. The singularities 1 — 4w 2 = 
also carry a x 1 ^ 2 behavior. Note that the order q of the linear differential equations of 
these $^''s follows a simple (parity dependent) linear rule as a function of n, namely 
q = (2n - (-1)" + 9) /4. 

Let us recall the "nickcllian" singularities (fT|) found by Nickel for the n— particle 
contributions \ to the magnetic susceptibility of Ising model pQ. These 
singularities, written as a function of Chebyshev polynomials of first kind 

h - Tfc ( cos (?)) +Tm ( cos (?))' °^ fc <™^ ( 42 ) 

can all be obtained from this "toy" model. Recalling the definition of the variable 
w = s/2/(l + s 2 ) which is self dual (i.e. invariant by s — > 1/s), one may introduce 

(n) 

the following partition on the singularities of the linear ODE corresponding to 3> y D . 
We denote by the polynomials where the roots (i.e. singularities) are given by 
(|42p. We denote by P^ the polynomials where the roots are not given by (|42|) but 
are in W c . Finally, we denote by -P 3 the polynomials where, at least, one root is not 
in W c - Along this partition, the polynomials for n = 3, • • • ,8, are respectively 

(void means there is no such type of polynomials): 

Pf ] = (1 - w)(l + 2w){l - l&w 2 ) ■ w, (43) 



for 



for 



P 2 (3) = void, P 3 (3) = l + 3io + 4w 2 



P 1 (4) = (1 - Aw 2 ) (1 - 16w 2 ) • w, (44) 
P 2 (4) = void, P 3 (4) = void 



PP = (1 + u,)(l - 16u; 2 )(l - 3w + w 2 )(l + 2w - Aw 2 ) ■ w, 

P 2 (5) = (1 -w)(l + 2w), (45) 

P 3 (5) = (1 - w - 3w 2 + 4w 3 )(l + 8w + 20w 2 + 15w 3 + Aw 4 ) 
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for 



for $g\ 



for $g\ 



pf 3 = (1 - w 2 ) (1 - 16w 2 ) (1 - 4w 2 ) (1 - 9w 2 ) ■ w, (46) 
f> = void, P 3 (e 



P, (6) = void, P 3 (6) = 1 - 10 w 2 +29w 4 



P[ 7) = w ■ (1 - 16w 2 )(l - 5w + 6w 2 - w 3 )(l + 2w- 8w 2 ~ 8w 3 ) 

(1 + 2w - w 2 - w 3 ), 
P 2 (7) = (1 + «j)(1 -3w + w 2 ){l + 2w- 4w 2 ), (47) 
P 3 (7) = (1 + 12u) + 54w 2 + 112w; 3 + 105w 4 + 35w 5 + 4w 6 ) 

(1 - 3w - I0w 2 + 35ui 3 + 5w 4 - 62w 5 + I7w 6 + 32w 7 - 16w 8 ) 



pf } = w ■ (1 - 16w 2 )(l - 4w 2 )(l - 2w 2 )(l - 8w 2 )(l +4w + 2w 2 ) 
(1 -4iu + 2w 2 ), 

P 2 (8) = (1- w 2 )(l-9w 2 ), (48) 
= l - 2Qw 2 + 2A2w A - 960u; 6 + 1685w 8 - 1138w 10 

for 

Discarding the singularities 1 — 16u> 2 = 0, one may remark, that for n odd, the 
singularities, roots of the polynomials P 2 , are also roots of the polynomials P[ n 2 \ 
The remark holds for n even, with the additional feature that 1 — Aw 2 = are common 
singularities to $>^\ n even. 



4-2. Singularities of the linear ODE's versus singularities of particular solutions 

Having not been able to find all the solutions of the linear differential equations 
satisfied by (|36[) . one can use the connection matrix method we introduced in [5] 
to find, for some values of n, the singular behavior of (|36|) at each regular singularity. 
This connection matrix method [8] was used to find the singular behavior of 
and X : the third and fourth "particle" contribution to the Ising magnetic 
susceptibility [S]. Let us sketch the method. The method consists in equating, at 
some matching points, the two sets of series corresponding, respectively, to expansions 
around a singular point and the nearest other singular point. The matching points 
have to be taken in the radius of convergence of both series. Connection matrices 
of two singularities that are not "neighbors" have to be deduced using some path of 
"neighboring" connection matrices. 

We denote by q the order of the linear differential equation. Connecting the local 
series-solutions at the regular singular points w = u>\, and w — u>2, amounts to 
finding the q x q matrix C(wi,W2) such that 

5 (lUl) = C{wi,w 2 )-S {w2) (49) 

where S^ 1 ^ and S^ 2 ^ denote the vectors whose entries are the local series-solutions at 
respectively the singular point ui\ and W2 ■ These solutions are evaluated numerically 
at q arbitrary points around a point belonging to both convergence disks of the scrics- 

(3) 

solutions. We have considered the linear differential equations corresponding to 3> D , 
$g } and 
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For the toy integral $> D whose linear differential equation has the singularities 
w = 0, w = 1/4, w = —1/4, w = 1, w — —1/2, w = oo and w = w\, u>2 where W\ 
and W2 are the roots of 1 + 3w + Aw 2 — 0, one obtains the following singular behavior 
around each singularity 



for the singularities 1 + 3w + Aw 2 = 0, namely b\ = 0. Recalling the singularity 
partition (|43p. one sees that only the singularities given by the polynomials where all 
the roots are in W c are actually present in the integral. 

We have also considered the linear ODE for which contains only singularities 
in W c . All these singularities occur in the integral. The next example is the linear 
ODE for All the singularities, given as roots of the polynomial p[ 5 \ occur in the 

integral, while the singularities roots of P 2 and P 3 do not. Note that two complex 
singularities, given by the polynomial 1 — w — 3w 2 + Aw 3 — 0, and all the roots of 
1 + 8w + 20w 2 + I5w 3 + Aw 4 = are not in W c . The last example we have considered 
is . Its linear differential equation has four singularities out of W c , namely the 
four roots of 1 — 10 w 2 + 29 w 4 = 0. Here also, we found that only the polynomials 
where all the roots are in W c , correspond to singularities of the integral. 



Note that the amplitudes at the singularities given by P2"]) can be obtained 

from 



where the sum is restricted to the integers that satisfy (|4"2"|) for the singularity w s . 
This formula is in agreement with the results of the connection matrix method, for 
n = 3, • • • , 6, confirming its efficiency. 

The natural question now, is to find out whether the Landau singularities 
identify, or have an overlap, with the singularities of the integral or, rather, with 
the singularities of the linear ODE. 

ff This is equivalent to the formula (14) forwarded by Nickel [T] for the x an d f° r the same type 
of singularities. 



$g } (singular, 1/4) = -— \n(x) ■ S ( 3 1/i] (x), 

47T 

(singular, -1/4) = -— -x ■ ln(x) • 4~ 1/4) (x), 

8 7T 

$g } (singular, 1) = -L ■ x' 1 ' 2 ■ S^\x), 
V6 

$g } (singular, -1/2) = 4= ' x ^' 2 ' s[~ 1,2 \x), 





(50) 
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4-3. Singularities from Landau conditions 

Before solving the Landau conditions for our integrals , some additional comments 
on Landau formulation should be underlined. The quantity Xj given in ([6]), is written 
aJ3 



^ - cosfa) + * J 1 - ( ±- - cos(& ) ) + / n ( 5 1 ) 



2 



Define 

mn(t\ = 

2w 



cos(C) = i-cos(^), (52) 



sin(C) = Jl-(J-- cosfoi)) + *0+ (53) 

and, similarly, cos(£„) and sin(£ n ), where <p is changed to (n — 1) </>. 
Using some identities on Chebyshev polynomials and noting that: 

z"" 1 ^) = T„_i (cos(C)) + i sin(C) • C/ n _2 (cos(C)) (54) 

the ratio p4|l . with the constraint (|35|) . is cast in the form 

cos(Cn) + T n _i (cos(C)) - i sin(C„) 4- i sin(C) ■ f/ n -2 (cos(C)) 



cos(Cn) - T n _i (cos(C)) - i sin(C„) - i sin(C) ■ C/ n -2 (cos(C)) 



(55) 



Some manipulations give for the integral 

$C = dio # ' cos(C n )-T n _ l( c s(C)) (56) 
where J7 n is the Chebyshev polynomial of second kind. In this form the denominator 
of the integrand is polynomial in w and cos(0), resulting from a "rationalization" 
of the integrand. This procedure introduces the "Galois companions" of the integral 
and we may expect the Landau conditions to generate, up to spurious singularities, 
all the singularities of the linear ODE rather than those of the integral. The second 
comment is on the branch points appearing in the numerator. The manifolds defining 
these singularities should be added in the Landau conditions. 

In the sequel, we give the singularities, obtained from the Landau conditions, 
in two cases. Firstly we allow the production of the singularities of the "Galois 
companions" , and consider the locus of singularities of the branch points. In a second 
step, we restrict the analysis on the integral without considering the branch points 
locus. 

The singularities, when the branch points appearing in the numerator are not 
considered in Landau conditions (see Appendix E), are the roots of: 



and 



TV, (i+l) +(-l)«-l-J- = (58) 
1 2w I 2w 



f i + is introduced for an easy control of our computations. 
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and the w-polynomials obtained by eliminating the variable v from: 

T„_ 1(W ) +T ^(±-v) -±= 0, (59) 

C/„_ 2 (v) - U n - 2 - vj = 

When the branch points are included, we obtain the additional polynomials (k and m 
are integers) , given by: 

^(i-^K-ir-^o, (so) 

and by the elimination of v from: 

T„-x(,)-(^±l)=0, ^-x(i-)±l = (61) 

Remarkably we find that the roots of the polynomials (jBTj) , (fS"5|) and (f5"5|) identify 
with the singularities of the differential equations governing (|36j) up to n = 8, i.e. 
with the singularities given in the polynomials P[ n \ P 2 an d • The check has 
been extented up to n = 16 with the singularities of the linear differential equations 
obtained, this time, modulo some primes. 

Note that the polynomials given by the first set (f57|) have all their roots in W c , 
and similarly for the second set (|58[) for n even. 

The polynomials (|60| , (|6Tj) give additional singularities not given by (|57ti59[) . For 
n = 3,4, 5, they read: 

n = 3, (1 + 3w)(l - 5w)(l - w - 4w 2 ) 

n = 4, (l±3«i-ro 2 )(l±6w+8w 2 ±4w 3 ) 

n = 5, (1 + 3 w) (4 w 3 + 3 w 2 - w - l) (5 w 2 + 5 w + l) 

(l - 16 w 2 - 2 + 56 w 4 - 16 w 5 - 63 w 6 + 8 w 7 + 16 w 8 ) 

(-l + 8w-20u; 2 + 17w 3 ). 

These singularities do not appear in the corresponding linear differential equation, 
they are spurious. 

Recall that these "spurious" singularities (with respect to the linear ODE) are 
obtained when the manifold of the branch points appearing in the numerator of (|56[) are 
considered in the Landau conditions. Let us now, exclude these (spurious) conditions 
on the branch points and, furthermore, let us control the solutions introduced by the 
normalization procedure in order to extract the " physical" Landau singularities^] from 
the "Galois" Landau singularities (see Appendix E). In this case, the singularities are 
obtained as roots of the gcd of the polynomials: 

l + E^J-L-ll =0 



and 



r^l^+l) +(-1)^-^ = 0, (63) 



1 By "physical", we mean "concerns" the integral, i.e. a particular solution of the linear ODE. 
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and the w-polynomials obtained by eliminating the variable v from 

T B _ l( tO -i-= 0, (64) 

U„-2 («) - Un-2 ( ^ - Vj = 0, 



I _ L / 1 \^ / I 



^fe-^j + v l -isr , v ■ u ^[^- v ) = ° 

We see that we obtain the same set of polynomials (|5T)) . (1551) and (|59|) dressed 
with constraints. The singularities, obtained from these conditions, correspond to the 
roots of the following polynomials 



n 


= 3, 


(1- 


4io)(l - 


f 2w)(l-w), 




n 


= 4, 


(1- 


4w)(l - 


f 4w)(l + 2iu)(l - 2w), 




n 


= 5, 


(1- 


4w)(l - 


f w)(l + 2w - 4w 2 ){l - 3w 


+ w 2 ), 


n 


- 6, 


(1- 


4w)(l - 


f 4w)(l - w 2 )(l - 9w 2 )(l - 


Aw 2 ), 


n 


= 7, 


(1- 


4w)(l - 


- 5w + 6w 2 - w 3 )(l + 2w- 


8w 2 - 8w : 






(1 + 2 


w — w 2 






n 


= 8, 


(1- 


4w)(l - 


f 4w)(l - 2w 2 )(l - 8w 2 )(l - 


-4w 2 )(l- 






(1 


- 4w + 


2w 2 )(l + 4w + 2w 2 ) 





9w 2 ) 

which are exactly (up to 1 + Aw for n odd) the roots given by P[ n \ i.e. Nickel's 
for the Ising model (|42|) . All these conditions (|62|l . ([63|l and (f64| give the "nickellian 
singularities" obtained from the simple form given in (|42p . 

Recall that we know exactly, by the results of the connection matrix method [8] 
given in the previous section, which singularities are actually present in the integrals 
for n — 3, 4, 5 and n — 6, and, especially, the non occurrence in the integrals of the 
singularities (not in W c ) roots of for n = 3, roots of P 2 ^\ P^ for n = 5 and 
roots of P^ for n — 6. So, and at least for n = 3,4, 5 and n = 6, the singularities 



obtained as roots of the polynomials (|62|) , (|63|) and (|64|) are in complete agreement 
with the singularities of the integrals. 

We thus see that, depending on how the Landau conditions are performed, the 
Landau singularities may be a set larger than the singularities of the linear ODE, or 
identify with the singularities of the ODE, or, even, could be the singularities of the 
integral. 



4-4- The singularities of the linear ODE's: an involved set with a structure 

Assuming that the polynomials ([57]) . ([55|) and (|55|) do indeed reproduce the 
singularities of the linear ODE corresponding to the integrals $^' ) for n > 16, one 
may use these formula to see how the singularities accumulate in the complex plane 
of the variable s, as n goes larger. Figure 1 shows the loci of the singularities outside 
the unit circle \s\ = 1. 
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Fig. 1 The roots of P 3 in the complex plane of the variable s, up to n — 47. 



These singularities with \s\ ^ 1 come only from the polynomials Ij58|) for n odd, 
and from the polynomials . Note that the singularities from {55J are all in the left 
half-plane of the variable s. The concentric circles are respectively (from the inner to 
the outer) |s| = 1/v2j l s l — 1 an d l s l = v2- The outer concentric circles correspond 
to the modulus in s of the quadratic numbers 1 + 3w + Aw 2 = (that read in the 
variable s, (1 + s + 2s 2 )(2 + s + s 2 ) = 0). We note the remarkable fact that all the 
singularities are lying in the annulus defined by these inner and outer circles and are 
far from the singularities s = ±1. These last two points seem to be accumulation 
points to which the singularities are converging as n goes higher. 

The proliferation of the singularities out of the unit circle \s\ = 1 is obviously 
polynomial. In fact, the degree of the polynomials Pg follows the simple rule 

deg(P 3 (ll) ) = ~ (2n 2 - 15) - { -^- (An + 1) , n = 3,4,--- 

When "all" the solutions of the ODE for each n are considered, their singularities 
form the unit circle \s\ = 1, and the structure around, shown in Fig. 1. 

4-5. More Kramers- Wannier duality subtleties: analytical continuation in s 

Let us consider, now, the linear differential equations corresponding to the integrals 
, which, written in the variable s, are obviously covariant by Kramers- Wannier 
duality s «-> \/s. Note that, this does not mean that the solutions of these self-dual 
differential equations are invariant by Kramers- Wannier duality s <-> 1/s. Consider 
first, the simplest integral of this family, namely < i ) ^' ) . As a function of the variable 
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w it is a solution of a third order differential operator which is the direct sum of two 
linear differential operators: 

Dv, © + ^ • Dw - ) (65) 

V (1 — lourjw; 1 — lour/ 

with general solution 

c + ci • iT(4w) + c 2 • CK{Aw) (66) 

where (CK) and if are respectively the (complementary) complete elliptic integral of 
first kind: 



CK(4:w) = K[y/l - 16 w 2 J (67) 

K(iw) = aFi ([1/2, 1/2], [1], 16u> 2 ) (68) 

The linear differential operator l]65p. written in the variable s, gives a third order 
linear differential operator L 3 (s) which reads: 

L 3 (s) = Ds © L 2 {s) where, 

L ^ - Ds2 + 8 \¥-iu • Ds - (irW (69) 

The third order linear differential operator L^{s) has the following general solution: 
c + ci ■ (1 + s 2 ) • if(s 2 ) + c 2 ■ (1 + s 2 ) • Cif (s 2 ) (70) 

The corresponding linear differential operators associated with the have the 

constant fonction as solution. Note that, for the integrals this constant^ can 

actually be calculated and is equal to 1/2. 
The particular solution $^ (s) reads: 

d>g } (s) = \ + I(l + s 2 )-if(s 2 ) (71) 
The general solution of £ 2 (s) can be written as: 

Mi-*g>(*)+ W^Q) (72) 

The two solutions of L 2 (s) are actually the same function for, respectively, the 
argument s and 1/s. Note that, this can be easily seen on the formal solutions 
of £ 2 (s) which take the same form 

*g>(u)-l/2, 

(*g> («) - 1/2) • m(u) + 1 (V + u 6 + |u 8 + • • •) (73) 

where u — s at the singular point s = and u = 1/s at the point at the infinity. 

The two quantities <&^(s) and <5^y(l/s) are two different solutions of the self- 
dual differential equations, in complete analogy with Y^(s) and Y^ 3 \l/s) of Section 
(3.3), especially their singularities which clearly break s <-» 1/s symmetry. 

We have here the exact equivalent of the situation we had [5] with the order-two 
linear differential equation corresponding to ■ The latter reads: 

£( 2 ) = 4w 4 ■ 2 Fi ([5/2,3/2], [3], 16u> 2 ) (74) 

(j We discard from now on, the factor 1/n! in the definition of the integrals 
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while the second independent solution \§\ reads in terms of the MeijerG function [28] : 



S 3 (w)= |MeijerG([[],[l/2,3/2]],[[2,0],[]],16^ 2 ) (75) 

Rewriting the order- two linear differential equation in the variable s, one obtains the 
general solutions: 

F(s) = d • i • ( 1 + 3 V 2 Fx ([5/2, 3/2], [2], i 



s^y v" ' " " " s 4 
+ c 2 - S 4 -(f + ,s 2 ). 2 ^([5/2,3/2],[2], S 4 ) 

- Cl -x (2) (l/s) +c 2 -x (2) (s) (76) 
The compatibility of (|74f and (|76|) corresponds to remarkable identities like: 

2*1(15/2, 3/2], [3],^i^) = (77) 
= (1 + U 2 ) 5 • ^([5/2,3/2], [2], U 4 ) 

or 

2 Fi ([1/2, 1/2], [1],^^) = (78) 

= (1 + U 2 ) • 2 F 1 ([1/2,1/2],[1], U 4 ) 

where u = s for small values of s and u = 1/s for large values of s. 

The previous calculations deal with simple linear ODE's and simple enough 
functions. However, for solutions of Fuchsian linear differential equations of higher 
order it becomes a wishful thinking to recognize, in the whole s-complex plane, that 
two functions are the "same function" for, respectively, the arguments s and 1/s. 

Consider, for instance, the integral $> D (s). The order four linear differential 
operator for & D (s) has the following direct sum decomposition: 

Ds®£ 3 (s) = Ds® (c 2 (s)-A(s)) (79) 

where the order-one, and order-two, linear differential operators are given in Appendix 
F. 

The integral $^(s) is the sum of the constant 1/2 and of a particular solution 
of£ 3 : 

S(C 3 ) = ci-5(A) +S(C 1 )- f ds^- (80) 



where S(Ci) is the solution of the self-dual order-one linear differential operator C\\ 
, . (l + s 2 ) (2 + s + 2s 2 ) , x 

S(Cl) = , , r- 81 

y/l + s + s 2 V2 - s + 2 s V2 + s + s 2 y/1 + s + 2 s 2 
The general solution of the order-two linear differential operator C 2 is a linear 
combination of Si(£ 2 ) and S 2 {L 2 ) given by: 

Si(£ 2 )(s) = -i- ' [Dk ■ K(s 2 ) + D E ■ E(s 2 )) , (82) 
S 2 (C 2 )(s) = 1 • S 1 (C 2 )(l/s) 
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where D , and De are polynomials in s given in Appendix F, and E(s 2 ) is the 
complete elliptic integral of second kind. The order-two linear differential operator £ 2 
is not self-dual, it is self-dual "up to s 2 " (see ([52]) ). 

To see that $^(s) and $^(l/s) are two different solutions of the linear 
differential operator (|79p . the integration in (|80|) has to be worked out. 

As far as the singular behavior of solutions of finite order linear differential 
operators with polynomial coefficients is concerned, the connection matrix method we 
introduced in [5] actually provides an answer to these difficult problems of analytical 
continuation in the whole s- complex plane. For the singular behavior of & D (s) at the 
points (2 + s + s 2 )(l + s + 2s 2 ) = 0, we have, typically, to compose three 3x3 connection 
matrices linking the point s = to si =—1/4+1 \fl / 4, root of 1 + s + 2s 2 = which is 
inside the unit circle, then to s = i, then to S2 = — 1/2 + iv7/2, root of s 2 + s + 2 = 
which is outside the unit circle. The corresponding singular behavior are: 

(singular, si) = b x ■ x~ 1/2 ■ S[ si) {x), 

$g } (singular, s 2 ) = b 2 ■ x~ 1/2 ■ sf 2) (x) (83) 

where S^ 3 ^ are series-solutions, analytical at x = (x is the expansion local variable 

(s ■) 

x = 1 — s/si), and such that S'j (0) = 1. The amplitudes at the singularities 

are b\ = and b 2 — (1 + i) ■ 7~ 1//4 /2. This means that the integral 4>^ is not 
singular at the two quadratic roots 1 + s + 2s 2 = which are inside the unit circle 
(|s| < 1). In contrast & D actually exhibits a singular behavior at the two quadratic 
roots 2 + s + s 2 = which are outside the unit circle (\s\ > 1). 

(5) 

We have also considered the singular behavior of Q D in the complex plane 
of the variable s. The singularities where at least one is out of the unit circle 
(|s| 1) are given by the roots of the polynomial which decomposes in 

two polynomials reading, respectively, (4 — 2s + 9s 2 — 2s 3 + 9s 4 — 2s 5 + 4s 6 ) and 
(4 + 6s + 7s 2 + 5s 3 + 2s 4 )(2 + 5s + 7s 2 + 6s 3 + 4s 4 ). The integral $^ is not singular at 
the roots of the first polynomial and at the roots of (2 + 5s + 7s 2 + 6s 3 + 4s 4 ) = 0, while 
its analytic continuation is actually singular at the roots of (4+6s + 7s 2 -|-5s 3 -|-2s 4 ) = 
which are exterior to the unit circle (|s| > 1). For Qp, the ODE carries as singularities 
out of the unit circle, the roots of 1 — 10w 2 + 29u> 4 = which reads in the variable s, 
16 + 24s 2 + 45s 4 + 24s 6 + 16s 8 = 0. The integral $^ is not singular at these points. 

As a result, one may imagine that the singularities of the $^(s), defined in the 
whole complex plane of the variable s by analytical continuation of high temperature 
s— series (valid for |s| < 1), actually correspond to the "nickellian singularities" (fTJ), 
together with some points in the "cloud" of singularities of Fig. 1 outside the unit 
circle (|s| > 1). 

In view of these results on the integrals and it is necessary to revisit, 

in the s variable, the connection matrix computation for the third contribution to the 
susceptibility of Ising model, given in [8]. 

We find that, similarly, to our toy integral X^ 3 ** is not singular at the two 

roots 1 + s + 2s 2 = 0, for which |s| < 1, while its analytical continuation is actually 
singular at the two roots 2 + s + s 2 — 0, for which |s| > 1. The corresponding singular 
behavior reads: 

X (3) (singular, si) = h ■ S[ Sl) (x) ■ x ■ ln(x), 
(singular, s 2 ) = b 2 ■ S^ix) ■ x ■ ln(x) 
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( s - ) 

where 5j are series-solutions, analytical at x = (x is the local variable x = 1- 



and such that S 1 , 
b 2 = (181V7 



(0) = 1. The amplitudes at the singularities are bi = and 
7i)/2567r. 

The three-particle contribution [8] to the susceptibility of Ising model, x^H 5 ); i s 
fundamentally a high-temperature quantity. The three-particle contribution, X^(s), 
together with x^0-/ s )i are different solutions of the self-dual linear differential 
operator for x^ 3 H s )- 

Similarly to the &^\s), one may imagine that the singularities of the x^ n H s )i f° r 
n odd, defined in the whole complex plane of the variable s by analytical continuation 
of high temperature s— series (valid for \s\ < 1), actually correspond to the "nickellian 
singularities" (p}, together with a similar "cloud" of singularities as in Fig. 1 outside 
the unit circle (\s\ > 1). 



4-6. Analytic continuation of Landau conditions in the variable s 

These Kramers- Wannier duality subtleties, and the breaking of the s <-*■ 1/s symmetry, 
can also be seen on the Landau conditions. Consider for instance, the polynomial 
1 + 3w + 4w 2 = which appears in the linear ODE for from the condition (|58|) 
and which does not occur in the integral due the constraint on this condition given in 
(ESJ). 

When one switches to the variable s, the polynomial (1 + 3u> 
(1 + s + 2s 2 )(2 + s + s 2 ). The constraint in (JgHJ) gives: 



Aw 2 ) becomes 



(i + s 2 ) (s - iy 



(84) 



2a - 



(1 + S 2 ) (S + I) 2 



(1 + s + s 2 ) 



The series expansions of the left-hand side of the algebraic condition 
the series expansions of: 



and 



(l + 2s) • (2 + s + s 2 ) = 



{s + 2) (2 s 2 + s + 1) = 



identify with 



(85) 



for respectively, small values of s and large values of s. The gcd in 
s, and for small values of s, then gives the two quadratic roots 2 + 



(86) 

written in 
s 2 = as 



("3") 

singularities of the integral $V> ■ These singularities are in the exterior of the unit 



circle (|s| > 1). Similarly, had the integral been defined for the large values of s, 



the gcd gives, as singularities, the roots of the polynomial 2 s 2 
in the interior \s\ < 1 of the unit circle. 



s + 1 = 0, which are 



5. Comments on the resolution of Landau conditions 

In Section (2), we have recalled some basic general ideas and definitions on the Landau 
conditions. We have used these conditions on our two families of integrals, for which, 
we knew, from the outset, the set of singularities occurring in the linear ODE and 
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in the integral. In practice the calculations can be slightly more subtle, and require 
cautiou jfjl . sometimes tricky and involved, analysis on one (resp. several) complex 
variable(s). For instance, the condition (f3"2")l in Section (|3.2p which rules out all the 
singularities not occurring in the integral, has to be worked out carefully to be able 
to control the ± signs floating around. 

The idea of the Landau conditions is to get some candidates for the singularities of 
some (multiple) integral from simple enough (algebraic) calculations on the integrand. 
In practice when the integrands are not rational expressions but algebraic ones, the 
algebraic character of the integrand introduces some ambiguity for each of the branchs 
cuts. 

For the calculations where the algebraic expressions are quite involved, the control 
of the Riemann sheet we stay on, may be quite tedious, sometimes hopeless. The 
algebraic integrand being solution of some polynomial, possibly of large degree, 
a possible approach amounts to considering all these roots together, namely the 
integrand of the actual "physical" (multiple) integral we are interested in, together 
with all his Galois "companions" . Recall that, as a consequence of the algebraic 
character of these integrands, the (multiple) integral we are analyzing is holonomic. It 
is a solution of a finite order linear differential equation with polynomial coefficients, 
and it is therefore, canonically associated with its "Galois companions" with respect 
to the linear ODE, namely the other solutions of the linear ODE. 

The calculations of the Landau conditions can now be performed for all these 
roots, leading to a larger set of Landau singularities that are probably good candidates 
for, not only the "physical" (multiple) integral we are interested in, but also for 
other mathematical expressions corresponding to (multiple) integral of his Galois 
"companions" integrands. 

The conditions, we have used to settle the set of singularities occurring in the 
integrals from the full set of singularities of the linear ODE, appear in (f3Tj) . (|32)l 
for the first family, and in (|103p . (|104p for the second family. Note that (|103|1 
is the denominator appearing in ()56|1 obtained with a "rationalization procedure", 
introducing the "Galois companions" of the integrand, and thus leading to the 
singularities of the linear ODE. We are aware that, for many examples (not of the 
Ising class we are interested in), the resolution of Landau conditions may be more 
tricky and/or more subtle. 

Let us close this section by a comment on the w versus s analytic continuation. 
We have used the self-dual variable w = s/2/(l + s 2 ) in our previous works (4j O El [8] 
on the susceptibility of Ising model since it appears naturally in the integrand ([3]) 
and, especially, it is much simpler as far as computational efforts are concerned. The 
differential equation for nas been obtained [4] with a series up to w 359 . In the 
variable s, a series up to s 699 would have been needed. When the distinction between 
the regimes of large values of s, and small values of s, is required, we have seen some 
of the subtleties that pop out. For the family, this subtlety appeared through a 
perfect square mechanism in the complex variable s. For the family, an equivalent 
mechanism of factorization appeared and the result depends on the regime considered 
(small s, large s). Note that this may not be the only mechanism that can be at work. 
One may imagine, for instance, a function, like the difference between the right-hand- 
side and the left-hand-side of ([77]). that evaluates differently according to the region 

ft See for instance the examples in Itzykson and Zuber's book |19j . or in Eden et. al [9]. 

U The Galois group of permutation of the roots of the polynomial is replaced here by the differential 

Galois group of the linear differential equation. 
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in the complex plane of the variable s we are considering. 



6. Towards Ising model integrals 

The two families of integrals presented in the previous sections are very rough 
approximations to the integrals (J5|). The first family considered the product on y i; 
integrated on the whole domain of integration of the fa. Here, we found a set of 
singularities occurring in the x an d the quadratic polynomial 1 + 3w + Aw 2 = 0. 
The second family was constructed, besides discarding the factor G^ n ' and the product 
on yi, by restricting the domain of integration on the principal diagonal of the angles 
fa. This resulted in a remarkable "memory" of the original problem since all the 
singularities ([42|) . forwarded by Nickel [HE], are reproduced, even the quadratic roots 
of 1 + 3w + Aw 2 = found (2 [5] for the linear ODE of x (3) - 

The following steps will be to incorporate, gradually, the various factors appearing 
in ([3]). One may, for instance, continue to discard the factor G^ n \ but incorporate the 
product on yj to the <S>p(w). 

In this move towards the n-particle contributions x^ °f the Ising model 
susceptibility, it is natural to consider the following family of integrals 

•s'-a-CnjC £)(n»)-f3fef 

which amounts to getting rid of the (G^) 2 in ([3|). 

One may restrict this multidimensional integral to a simple integral on one angular 
variable in the same way we did it for our previous "diagonal toy model" in Section 
©: 

n\ J Q 2?r 

For this last family, one obtains for n = 3, the singularities found for the linear 
ODE of the diagonal model integral of Section and also the singularities 

(1 + 3u?)(l — 5w) = 0. These singularities emerge as a consequence of inclusion 
of the yi. Recall that these last singularities do not occur [3] in the linear ODE 
corresponding to x^- Do they disappear in x^ as a consequence of the enlargment 
of multidimensional integration domain, or as a consequence of the inclusion of the 
G(") factor ? 

Note that, as we move closer to the original integrand and integral ©, for 
instance considering the multidimensional integral (|57| . the calculations 4, 5, 6J (series 
expansions, search of the linear ODE) become much harder. Nevertheless, an exact 
knowledge of the singularities of the linear ODE's can be achieved using a brand new 
strategy that amounts to getting very large series for these integrals modulo primes, 
and in a second step, get the corresponding linear ODE's also modulo prime^. Ideally, 

K This can be seen as a new method for obtaining exact results in lattice statistical mechanics, and, 
in any case, extremely large series expansions. In forthcoming publications with A. J. Guttmann, 
I.Jensen and B. Nickel we will detail such "extreme" calculations of lattice statistical mechanics 
modulo the primes providing more than 1600 coefficients of the high and low temperature series 
expansions of the susceptibility \ of the square Ising model, as well as more than 2000 coefficients 
for x' 5 ', and, for given primes, with 6000 coefficients. 
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for a large enough set of such calculations modulo different primes, one can get (from 
a Chinese remainder procedure) the exact linear ODE we are seeking for. In the case 
where the number of primes and calculations modulo these primes are not numerous 
enough to build the linear ODE, they are sufficient to be sure that a certain polynomial 
(like the ones we displayed in Section actually occurs in the head polynomial of 
the linear ODE, or discard its occurrence. 

7. Conclusion 

We have designed two families of "Ising class integrals" (Y^(w), 3>^\w)), for which 
we obtained the Fuchsian linear differential equations for large enough values of n. 
Each Fuchsian linear differential equation provides a set of singularities. We have used 
a direct resolution (for the first family of integrals) and our connection matrix method 
[5] (for the second family of integrals) to obtain the subset of singularities occurring in 
the integrals. For both families of integrals, we solved the Landau conditions and found 
that the singularities we obtained may identify, or even extend, the singularities of the 
linear ODE or the integrals. While it is obvious that the linear ODE, which has this 
holonomic integral as solution, gives all other solutions with their singular behavior, 
it is remarkable, as far as the locus of singularities is concerned, that the Landau 
conditions calculations on the integrand give informations on the other singularities 
not carried by the integral. 

For these "Ising class integrals" , we found that each family of integrals is singular 
in the domain W c . Our integrals are defined for small values of the variable w and 
elsewhere by analytical continuation. This variable w = s/2/(l + s 2 ) behaves in 
an equal footing for small and large values of the variable s. Switching to this last 
s variable, we showed that the analytical continuations (from small values of s for 
which the integrals are defined) of these two families of integrals are singular on points 
occurring on the unit circle \s\ = 1, but also on some points in the exterior of this 
circle (|s| > 1). 

In view of these results on our toy integrals, we revisited the singular behavior 
of the third contribution to the magnetic susceptibility x of the Ising model. We 
also found here that i defined as series expansions for small values of the variable 
s and analytically continued to large values of s, is actually singular outside the unit 
circle, |s| > 1, at the points 2 + s + s 2 = 0. If a similar scenario remains valid for the 
other X^( s )i f° r n odd, one may imagine that the singularities of the X^ n \ s )i f° r n 
odd, defined in the whole complex plane of the variable s by analytical continuation 
of high temperature s— series (valid for |s| < 1), actually correspond to the "nickellian 
singularities" ([1]), together with a similar "cloud" of singularities as in Fig. 1 outside 
the unit circle (\s\ > 1). 

Our two families of integrals are introduced to mimic the n— particle contribution 
X~ n ' to the susceptibility of the square Ising model. In a forthcoming publication we 
will see to what extend these two families of integrals, and similar toy models, give 
information on the singular behavior of these integrals in the complex plane of the 
variable s and on the singularities that might occur in the X ■ 
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8. Appendix A 

We give in this Appendix the order six differential equation corresponding to 

= (89) 



that reads (F(w) = Y^(w)) 

6 ,„ 

£ ««h • ^r F H = ( 9 °) 



where 



with 



n=0 



a 6 (w) = (l-5w)(l + 3u;)(l + 3w + 4w 2 )(l + 4w) 2 (l-4w7) 4 

(l-io) w 6 -P 6 (w) 
a 5 (w) = 2 (l + 4w) (-l + 4w) 3 w 5 ■ P 5 (w) 

a 4 {w) = 2(1 -Awfw^ ■ P 4 (w), a 3 {w) = -4 (1 - 4 w) ■ w 3 ■ P 3 {w) 
02 (tw) = 4 w 2 ■ P2(w), ai(w) = 8w ■ Pi (to), ao(w) — P (w) 

P 6 (w) = 464889446400w 17 - 5967097036800 w 16 + 2518865879040 w 15 
+ 1205147090944 w 14 - 209927829504 w 13 - 23192809472 w 12 

- 94576000000 w 11 + 15333496352 w 10 + 6474551520 w 9 + 434664056 w s 

- 768760124 w 7 + 47392381 w 6 + 24032176 w 5 - 3740777 w 4 
+ 106680 w 3 + 7443 w 2 - 364 w + 49 

P 5 (w) = 5020806021 120000 w 24 - 67498529469235200 w 23 
+ 30557536675430400 w 22 + 38088998121308160 w 21 

- 4237195746729984 w 20 - 9069732118712320 w 19 - 2846232832169984 w ls 
+ 241117241392640 w 17 + 618725352473088 w w + 201166769452896 w 15 

- 61245521268224 w 14 - 22970672645976 w 13 + 1946216478788 w 12 
+ 1256316716689 w 11 - 2295775019 w w - 38685730792 w 9 

+ 7182763949 w 8 - 2662652062 w 7 + 133819572 w 6 + 107083800 w 5 

- 16853044 w 4 + 499141 w 3 + 34119 w 2 - 784 w + 147 

Pi(w) = 562330274365440000 w 2e - 7836264826719436800 w 25 
+ 3202063539476889600 w 24 + 3620634209723351040 w 23 

- 548610793223356416 w 22 - 626862114022522880 w 21 

- 294517075356614656 w 20 + 30188250588135424 w 19 
+ 71955514644985856 w 18 + 6645859876651520 w 17 

- 7980786065196160 w 16 - 956452174732544 w 15 
+ 546200976422432 w 14 + 127674138396992 w 13 

- 45199723462020 w 12 - 9028528613608 w 11 + 2584065207968 w w 
+ 428924564347 w 9 - 59295462952 w s - 29452365937 w 7 

+ 3600235371 w 6 + 631583771 w 5 - 118172923 w 4 + 3790389 w 3 
+ 186393 w 2 - 5138 w + 931 
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P 3 (w) = 2811651371827200000 w 27 - 41266384607379456000 w 26 
+ 24773736838201344000 w 25 + 1 1079955148990054400 w 24 

- 4824757963497799680 w 23 - 848931834754990080 u; 22 

- 1254549778859409408 w 21 + 332436165839831040 w 20 
+ 230844285233797120 w 19 + 23397830675772928 w 18 

- 35731466794253184 u; 17 - 7015978530910848 w 16 
+ 3334508426832544 w; 15 + 129769752723552 w 14 

+ 64823043832572 u; 13 + 9687172326860 w 12 - 27269045355036 w 11 
+ 668812543530 w; 10 + 1232145138001 w 9 + 217652915670 w s 

- 113666172030 w 7 + 2636314376 w 6 + 2843958768 w 5 

- 343848072 w 4 + 7774186 w 3 + 533670 w 2 - 21133 w + 2254 
P 2 (w) = 84349541 15481600000 w 28 - 130054335243485184000 w 27 
+ 101906251242799104000 w 26 + 14913046820133273600 w 25 

- 16311072130508390400 w 24 + 3172467314793676800 w 23 

- 7523898588462022656 w 22 + 4342034847529598976 w 21 
+ 349567990475905024 w 20 - 1019978058716864000 w 19 
+ 135175191464703616 w 18 + 154325934113136384 w 17 

- 39182557846936736 w 16 - 6064402640954016 w 15 

- 721601939503388 w 14 + 778716878219200 w 13 + 414373001033428 u 

- 145383247830248 w 11 - 7935890967120 w 10 + 857310545482 w 9 
+ 2441550352242 w & - 452151200437 w 7 - 32488744995 w 6 

+ 15597930449 w 5 - 1389870493 w 4 + 21276545 w 3 
+ 2254771 w 2 - 106015 w + 7987 

Pi (to) = 421747705774080000 w 28 - 6710038906798080000 w 27 
+ 4873843971588096000 w 26 + 1036215739573862400 w 25 

- 616703179652628480 w 24 ~ 120979270258114560 w 23 
+ 364788984254619648 w 22 - 1023471173198333952 w 21 

- 201098747980344320 w 20 + 953808829975566464 w 19 

- 246857972154239232 w ls - 152394614023808544 w 17 
+ 55269675209779600 w 16 + 3979478375375692 w 15 

+ 612261498482876 w 14 - 759726033790536 w 13 - 527516229973280 u 
+ 188404041710426 w 11 + 5343977363342 w 10 - 599703686985 w 9 

- 2917831385349 w 8 + 538227200312 w 7 + 39103736504 w 6 

- 18721639514 w 5 + 1671294354 w 4 - 25882894 w 3 

- 2676882 w 2 + 122171 w - 9261 

PoH = -32133158535168000 w 26 + 695978420915404800 w 25 

- 544662214410240000 w 24 + 3609286260943749120 w 23 

- 7536908069936824320 w 22 + 10191253169267146752 w 21 
+ 2316960321752137728 w w - 8808777472942522368 w 19 
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+ 2190279912630165504 w 18 + 1472339197313387520 w 17 

- 535796399651893248 w 16 - 40698424759766784 w 15 

- 1761836408896896 w 14 + 6557432066372256 w 13 

+ 5053923219182592 w 12 - 1884372068029344 w 11 - 14699422810656 w 10 
+ 1982678836896 w 9 + 27382956586968 w 8 - 5031984639480 w 7 

- 368590752120 w 6 + 176119304760 w 5 - 15746069112 w 4 
+ 246121272 w 3 + 24962808 w 2 - 1111320 w + 84672 

9. Appendix B 

In this Appendix, we derive the yW expansion ([!§]) , At the first step, the n — 1 
dimensional integral is translated to a n dimensional one, the angular constraint ([7]) 
being taken into account by a Dirac delta distribution 8(527=1^) wrnc h is Fourier 
expanded as (with Z p = exp(i</> p )) 

n oo . 

2tt ■ 8(J2 <Pi) = ( Z * Z * " ' Z «) ( 91 ) 

i— 1 k— — oo 

The integrals (fT4|) read 



00 , /-2ir 7i 



y (92) 



fc=— 00 

From the Fourier expansion (|16[) , one has the trivial integration rule 



^g(2/ 2 ^ fc ) = ±w 2 -D{k) (93) 







One gets 

00 

r (n) H = (4w 2 ) n - D ( k ) n ( 94 ) 

fc— — 00 

and with the definition ([T7]), we finally obtain p^|) . 
10. Appendix C 

In this Appendix we give for the y 2 — product model, the differential equation for 
n = 4 and the polynomials Ri and R 2 corresponding to the solution Y^(w). The 
differential equation for (w) reads 

a 2 ^Y^ + ai ^Y^ +a Y^ = (95) 
aw z aw 

where 

a 2 = (1 - 16 w 2 ) 2 (l-4w 2 ) (l-20w 2 + 16w 4 -16w 6 ) w 2 ■ P 2 (w), 
a x = 2w (l - 16w 2 ) Pi(w), ao = 8P (w) 
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P 2 = -15 + 1086 w 2 - 28106 w 4 + 328716 w 6 - 1946216 w 8 
+ 4791440 w 10 - 4697088 w 12 + 8682368 w 14 - 17308416 w 16 



29 



1781760 w 



18 

Pi = 75 - 7716 w 2 + 319862 w 4 ~ 7044848 w 6 + 92471304 w " 

- 761556288 w 10 + 3861020800 w 12 - 11558422912 w 14 
+ 2206951 1424 w w - 42440728576 w 18 + 101426089984 

- 147848548352 w 22 + 72672886784 w 24 - 17528389632 
+ 2736783360 w 28 , 
P Q = -45 + 5400 w 2 - 261002 w 4 + 6787240 w 6 - 107032250 w 8 
+ 1079592480 w 10 - 6995904504 w 12 + 28478625024 w 14 

- 72563446848 w 16 + 142196898048 w 18 - 402528466944 w 20 

+ 1067417646080 w 22 - 1296486492160 w 24 + 418540593152 w 26 

- 55343579136 w 28 + 136 8 39 1 68 00 w 30 . 

For the solution Y^(w), the polynomial P^ and P^ are given in (|23|) and 

r[ 5) = 15 - 950 w 2 - 550 w 3 + 25581 w 4 + 14856 w 5 - 363086 w 6 
- 57174 w 7 + 2727087 w 8 - 1243960 w 9 - 9109540 w 10 
+ 10083482 w 11 + 8271200 w 12 - 23345456 w 13 - 2557976 w 14 
+ 19604672 w 15 + 3295040 w 16 - 1531392 w 17 - 513024 w 18 , 

R ( 2 5) = 3 + 6 w - 103 w 2 - 292 w 3 + 1207 w 4 + 2214 w 5 - 7360 w 6 
+ 3674 w 7 + 6896 w 8 - 21328 w 9 - 13080 w 10 - 96 w 11 - 256 w 12 . 

11. Appendix D 

(3) 

The differential equation satisfied by § D reads: 

n.= 1 



l W 20 
>w 26 



where 
with 



n=l 

= (1 - w) (1 + 2w) (1 - 4w) (1 + 4w) (1 + 3u> + 4w 2 ) W 2 • P±(w), 
a 3 (w) = w-P 3 (w), a 2 (w) = P 2 (w), a x (w) = P x (w) 

Pi{w) = -18 - 126 w - 536 w 2 - 581 w 3 + 11332 w 4 + 56216 w 5 
+ 141103 w 6 + 316146 w 7 + 324516 w 8 - 102512 w 9 + 512 w 10 
+ 104448 w 11 , 

P 3 (w) = -36 - 450 w - 270 w 2 + 17199 w 3 + 110892 u; 4 + 331122 w 5 

- 365893 w 6 - 9173304 w 7 - 40917443 w 8 - 92069955 w 9 

- 128675122 w 10 - 89628548 w 11 + 223226064 w 12 
+ 725436224 w 13 + 509586688 w 14 - 185729024 w 15 
+ 25788416 w 16 + 133693440 w 17 , 



class 30 



w 5 
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P 2 {w) = 36 + 396 w + 61 14 w 2 + 59550 w 3 + 302337 w 4 + 857367 
+ 887406 w 6 - 6754260 w 7 - 39314109 w 8 - 48807357 w 9 
+ 175442484 w 10 + 683394108 w 11 + 1327431600 w 12 
+ 2198750784 w 13 + 1607543040 w 14 - 479219712 w 15 
+ 35880960 w 16 + 320864256 u; 17 , 

Pi(w) = 108+ 1512 w + 24282 w 2 + 230094 w 3 + 1072500 
+ 3034020 w 5 + 6713730 w 6 + 7870578 w 7 + 26737134 ^ 
+ 253594938 w 9 + 799267644 w 10 + 1253301264 w 11 
+ 1490308224 w 12 + 1015581696 w 13 - 338632704 w 14 
- 1179648 w 15 + 160432128 w 16 . 
The differential equation satisfied by $^ reads: 

71—1 

a 4 (w) = w 3 {l + 4w)(l + 2w)(l-2w){l-4:w)-P 4 (w), 
a 3 (w) = 2w 4 -P 3 (w), a 2 (w) = 3w ■ P 2 (w) , a x (w) = 3 • P x (w) 

P 4 (w) = l + 19w 2 + 314w 4 + 512u; 6 , 

P 3 (w) = -81 - 1010 w 2 - 1756 w 4 + 89408 w 6 + 114688 w s , 

P2(w) = -1 - 53 w 2 + 98304 w 10 + 117888 w s + 742 w 4 + 24896 w 6 , 

P 1 (w) = 1 + 53 w 2 + 2042 w 4 + 15616 w 6 . 

12. Appendix E 

In this Appendix, we give the proof for the solutions of the Landau conditions. One 
may use the form of the integrand given in (|56|) , we instead use another form to show 
clearly how the singularities of the "Galois" companions are reproduced. 
The variable Xj given in ^ from the definition (|52[) can be written as: 

Xj{4>j) = exp(iCj) (98) 
The integral (f3"6")) then becomes 

D n\ J 2ir tan(V>/2)' 
for which, the Landau conditions read: 



where 
with 



- - I oZ nrTrrToT' ^=(n-i)c+& (99) 



ai ■ ip = 0, a 2 ■ 4> — 0, 

a 3 -sin(C) 2 = 0, a 4 -sin(C„) 2 = 0, 

ai-i>' + a 2 + 2a 3 ■ sin(^) cos(C) (100) 
+ 2ai4 • sin((n — 1)0) cos(( n ) = 

where the parameters ctj, i = 1, •••4, should not be, all, equal to zero. In the sequel, 
we give the proofs along the possible values (zero or not) of the two parameters a 3 
and «4. 
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12.1. a 3 = a A = 

In this configuration, we are discarding the conditions on the branch points. The 
Landau conditions for ct\ and cti not both equal to zero give: 

C„ + (n - 1)C = 0, 0=0 (101) 

or: 

Cn + (n-l)C=0, sin((n ; 1)0) = sin^ 

sm(Cn) sm(C) 

These equations correspond (with sin(£) ^ 0), respectively, to an end-point singularity 
and to a pinch singularity. The condition QlOip being contained in (| 102[) for = 0, 
one may consider the last. 

The condition + (n — 1)( = is equivalent to 

cos(C„) = cos ((n -1)0, (103) 
sin(C„) = -sin((n-l)C). (104) 
Using |52|) for cos(( n ) written as 

cos(C„) = ^- - cos((n - 1)0) = -L - T„_! (cos(0)) (105) 
one obtains for (|103[) : 

i-T n _i(cos(0)) = cos((n-l)cos- 1 ( -^-cos(0))) (106) 

Denoting v = cos(0), and by definition of the Chebyshev polynomials, one gets 

i-T^(v) = T^(±-v) (107) 
which is the first condition given in (|64p. One has 



sin ||;/ - DO = \/ | - ( -L - CO s(0) j • C7 n _ 2 - COS(0) 



and taking sin(^„) as 



sin(C„) = y/1 - cos(C„) 2 = ^1- (^-T„_!(cos(0))) 2 (108) 
the condition (|104p gives 

l-(^-r„-i(«)) 2 (109) 



which is the third condition given in (|64p and is a constraint to (JTUTJ). 

The second condition in (| 1 02|) is automatically satisfied for = and = n. 
Plugging these values in (|107j) and (|109p gives respectively 

1 +c, »- 2 (i- 1 ) = °> < 110) 
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which are the constraints given in (|62|) and (|63|) . The last have solutions only for n 
even. 

At this point, the second condition in (|102p has not been used. Using (1104p in 
this last condition, one obtains 

sin((n- 1) (cos- 1 ^)) sin (cos" 1 ^)) = 

sin((n- l)cos- 1 (l/2w- z)) sin (cos- x (l/2u) - «)) 1 J 

yielding: 

C/„- 2 (u) -^-af^-^) = ( 112 ) 

which is the second of equation in (|64p . 

Let us sum up for this case. Using the first condition in (|102p . giving both 
conditions (|103|) and (|104p . leads to a smaller set of singularities. While imposing 
the only condition (|103[) allows the appearance of the singularities of the "Galois 
companions" . 

12.2. a 3 ^ 

In this configuration, if = 0, the condition sin (£) = gives directly w = ±1/4. 

For «4 7^ 0, we have sin (C) = sin(£„) = 0, which, using the definition of cos(£) 
and cos (£„) leads to 

r^(-L -(-!)*) + (-!)« --L = „, (113) 
where k and m are integers. 

12.3. a 3 = and a 4 ^ 

In this configuration, additional polynomials are obtained for ct\ ^ resulting in 

(n-l)C +Cn = 0, with, cos(C) 2 = 1 (114) 

leading to 

T n _ l(u) -(i- ± l) = 0, and ^(-L-t,) ±1=0 
which is solved by elimination of v. 

13. Appendix F 

The order one and order two differential operators occurring in the order-four linear 
differential operator 

Ds@C 3 (s) = Ds® (c 2 (s) -£i(a)) (115) 

(3) 

written in the variable s and corresponding to $> D read: 

= Ds-\QL, C 2 (s) = Ds* +^.Ds+^ 
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with 

Qi = (1 + s 2 ) (s 2 + s + l) 

x (2 s 2 - s + 2) (2 s 2 + s + 1) (s 2 + s + 2) (2 s 2 + s + 2) , 
Qo = (s 2 - 1) 

x (8 s 8 + 12 s 7 + 26 s 6 + 19 s 5 + 32 s 4 + 19 s 3 + 26 s 2 + 12 s + 8) 



and 



.21 



P-2 = s (1 + s 2 ) {2s 2 + s + 2) 2 (s 2 - l) 2 (2s 2 + ,s + 1) (,s 2 + s + 2) , 

(s 2 + s + 1) (2 s 2 - s + 2) (64 s 16 + 328 s 14 + 596 s 13 + 4434 s 12 

+ 6621 s 11 + 15377 s w + 16897 s 9 + 22822 s 8 + 16897 s 7 + 15377 ; 
+ 6621 s 5 + 4434 s 4 + 596 s 3 + 328 s 2 + 64) , 

Pi = (s 2 - 1) (2 s 2 + s + 2) (l + s 2 ) (2048 s 28 + 38 40 s 27 + 177 92 s 26 

+ 63552 s 25 + 38 1 504 s 24 + 12 1 78 88 s 23 + 33 5 1 080 s 22 + 6980836 ; 
+ 12444792 s 20 + 18 5 5 6 93 s 19 + 23711235 s 18 + 26152837 s 17 
+ 24357824 s 16 + 19695118 s 15 + 12653972 s 14 + 7020790 s 13 
+ 2315882 s 12 + 789433 s 11 - 161991 s 10 + 312961 s 9 + 197622 s 8 
+ 343348 s 7 + 123248 s 6 + 49472 s 5 - 28512 s 4 - 28608 s 3 
-17152 s 2 - 3840 s - 1024), 

P = 2048 s 33 + 3072 s 32 + 17152 s 31 + 69760 s 30 + 526400 s 29 
+ 1787232 s 28 + 50 9 4 2 88 s 27 + 11072344 s 26 + 19 78 1 548 s 25 
+ 27467354 s 24 + 26 98 1 054 s 23 + 5 1 23 6 35 s 22 - 46141730 s 21 

- 134419668 s 20 - 24 2 2 3 9 80 s 19 - 354116927 s 18 - 430640772 s 17 

- 463578366 s 16 - 433931810 s 15 - 371127411 s 14 - 281456218 s 13 

- 203391640 s 12 - 133972232 s 11 - 89380041 s 10 

- 55343916 s 9 - 34593136 s 8 - 19033912 s 7 - 9514544 s 6 - 3824384 s 5 

- 1169344 s 4 - 213760 s 3 - 768 s 2 + 13824 s + 2048. 

The polynomials Do, Dk and occurring in the solutions (f82|) of the order 
two linear differential operator £2(3) are: 

Do = s (s - 1) (1 + s) (1 + s 2 ) (2 + s + 2 s 2 ) (2 - s + 2 s 2 ) 

(1 + s + s 2 ) (2 + s + s 2 ) (1 + s + 2 s 2 ) , 
£» x - (s-l)(l + s) 2 (l + s 2 ) 2 

(8 s 8 + 4 s 7 + 22 s 6 + 9 s 5 + 93 s 4 + 60 s 3 + 100 s 2 + 32 s + 32) , 
D E = 2 (1 + s 2 ) (1 + s) 2 (2 + s + 2 s 2 ) 

(8s 6 - 2s 5 + 17s 4 - 10s 3 + 17s 2 - 2 s + 8) . 
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